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Abstract. We consider dynamical decoupling schemes in which the qubit is 
continuously manipulated by a control field at all times. Building on the theory 
of the Uhrig Dynamical Decoupling sequence (UDD) and its connections to 
Chebyshev polynomials, we derive a method of always-on control by expressing the 
UDD control field as a Fourier series. We then truncate this series and numerically 
optimize the series coefficients for decoupling, constructing the CAFE (Chebyshev 
and Fourier Expansion) sequence. This approach generates a bounded, continuous 
control field. We simulate the decoupling effectiveness of our sequence vs. a 
continuous version of UDD for a qubit coupled to fully-quantum and semi-classical 
dephasing baths and find comparable performance. We derive filter functions for 
continuous-control decoupling sequences, and we assess how robust such sequences 
are to noise on control fields. The methods we employ provide a variety of tools 
to analyze continuous-control dynamical decoupling sequences. 



1. Introduction 

Quantum information processing seeks to exploit quantum mechanical systems to 
store and manipulate information in novel ways. Typically, however, this quantum 
behavior exists at the microscopic level, such as a single atom or electron, making 
control by the experimenter challenging. Moreover, such a system is never perfectly 
isolated. Stray electromagnetic fields, thermal fluctuations, etc. can perturb the 
quantum information carrier, imparting noise into the system. This noise can disrupt 
the quantum information processor and negate any benefits to operating in a quantum 
system. 

Quantum states can be protected from noise with active control fields. Before 
the field of quantum information emerged, a wide variety of techniques for dynamical 
decoupling (DD) of a spin ensemble were developed in the nuclear magnetic resonance 
(NMR) community. The first discovery was the Hahn spin echo sequence [l]. By 
repeating the spin echo sequence several times in succession, one finds the Carr- 
Purcell (CP) sequence (2). This sequence was soon improved by selectively choosing 
the axes around which one rotates the spin vector, yielding the workhorse Carr- 
Purcell-Meiboom-Gill (CPMG) sequence J3j. Ever more complicated sequences were 
developed in ensuing decades, allowing detailed characterization or manipulation of the 
environments of nuclear spin ensembles, or the compensation against inhomogeneous 
control fields. 
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Dynamical decoupling found a renewed purpose in quantum information science. 
For example, quantum error correction systems [4 10 demand an error-per-operation 
of less than about 1% [7UTT] . Such low error rates can be difficult to achieve in practice 



since many candidate physical systems will lose their qubit coherence very quickly 
owing to environment-induced decoherence; as a result, quantum error correction may 
fail due to short qubit lifetimes. Dynamical decoupling schemes offer a relatively 
simple solution to this problem: by applying a determined sequence of control pulses, 
one can significantly enhance the lifetime of a qubit. 

Early work in dynamical decoupling sought to suppress decoherence by applying 
a periodic sequence of instantaneous "bang-bang" pulses, which periodically flip the 
state of a qubit and undo the coupling to the bath 12 15 . This work was extended to 
Hamiltonian engineering by using decoupling sequences to selectively enable coupling 
Hamiltonians 
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Subsequently, attention turned to using Eulerian graphs to 
design sequences using bounded control operations that were robust to many types 
of systematic control errors 



gates 
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This led to the notion of dynamically-corrected 
which combine techniques from dynamical decoupling and composite pulse 
sequences from NMR to produce error-suppressing quantum gates. 

The idealized "bang-bang" control pulse is instantaneous, like the Dirac delta 
function. However, real physical pulses can only approximate "bang-bang" control, 
because such idealized control would require infinite power. The implications of a real, 
continuous-time function can be significant. Most dynamical decoupling sequences 
are designed to correct noise errors between pulses, but offer no intrinsic protection 
to errors during pulses. Early studies of always-on dynamical decoupling arose in 

More recently, optimal 
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the problem of dipolar decoupling in solid-state NMR 
control approaches to customizing continuous, bounded controls to decouple a specific 
but arbitrary noise bath were considered in Refs. 21 22 . Of particular relevance to 



the present work are approaches that extend optimized bang-bang sequences such as 
the Uhrig Dynamical Decoupling (UDD) sequence [23], which is derived by canceling 
general noise Hamiltonians order-by-order in perturbation theory. Extensions of UDD 
to continuous and bounded control were studied by Uhrig and Pasini 24|25 . This work 
introduced pulses that correct errors during their own evolution to finite order in a 
time-dependent perturbation expansion. These sequences are termed "realistic UDD," 
or RUDD. Subsequent work showed how to produce arbitrary gates decoupled from 
environment noise using bounded-strength controls [26j|27], as well as the application 
of modifications based on optimal control given known system drifts 28 . Dynamical 
decoupling with such "finite-width" pulses represents efforts to improve existing 
sequences by accounting for the real continuous-time nature of control fields. 

Dynamical decoupling sequences where the control field is always or almost 
always on pose some unique challenges for sequence engineering and characterization. 
Methods for characterization that are exact in the "bang-bang" limit of DD and 
approximately correct with short but finite pulse widths require continuous analogs. 
In this work, we examine methods to both create and characterize always-on sequences, 
by considering continuous extensions of the UDD sequence. We introduce a new DD 
sequence called CAFE (Chebyshev and Fourier Expansion), and compare it to similar 
RUDD sequences. 

The paper is organized as follows. Section [2] gives a derivation of the theory 
behind CAFE, presenting it as a continuous extension of the bang-bang UDD 
sequence. Ultimately, CAFE is derived from a particular heuristic rather than from a 
fundamental determination of optimality, but this is typical of continuous sequences. 
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Critical for evaluating the heuristic, then, are methods of characterizing the sequence, 
and careful choices of sequences for comparison, as discussed in section [3] Section [4] 
tests the CAFE sequence as well as RUDD sequences for performance in decoupling 
a qubit from a simulated quantum bath, represented by a small dipole-coupled spin 
lattice. Section [5] studies the frequency response of CAFE and compares analytical and 
simulated filter functions for this sequence and RUDD sequences. Section [6] simulates 
the scenario where a qubit is coupled to semiclassical noise and where control fields 
for decoupling are noisy, again for both CAFE and RUDD sequences. 

2. Derivation of CAFE sequences 

We consider a qubit that is exposed to an environment that causes phase errors, or a 
dephasing bath. For simplicity, we parameterize our system in terms of Pauli operators 
a 3 , so that the phase errors correspond to o~ z operators and the control of the qubit 
is proportional to the o~ x operator. The total system Hamiltonian is 

H(t) = 1 ® B + a z <g) B z + a{t)a x ® 1, (1) 

where 1 refers to identity in its subspace. In this expression, 1 ® Bq is a pure bath 
evolution term, a z ® Bz is the qubit-bath coupling term, and a(t)a x ® 1 is the time- 
dependent term representing the experimenter's control. In the following analysis, we 
find it useful to define the time-integrated control function 

/3(f) = 2 [ a(s)ds. (2) 
Jo 

The quantity /3(f) has an intuitive meaning — it represents the total angle (commonly 
called "pulse area") that the qubit state is rotated on the Bloch sphere, in this case 
around the X axis. We now transform our system to the interaction picture, so that 
the interaction Hamiltonian becomes 

n iut (t) = cos \p(t)] <r z <g> B(t) + sin [/3(f)] a v ® 5(f), (3) 

where 

B(t) = e lBat B z e- lBot . (4) 

We see that, in the interaction picture, the system-bath coupling terms are 
functions of /3(f). In the case of "bang-bang" 7r-pulse control, /3(f) is a piecewise- 
constant function that is everywhere a multiple of n, so that sin [/3(f)] = everywhere 
and we can replace cos [/3(f)] with a "switching function" y(f) that takes on the values 
±1. For example, this is the formalism considered in Ref. [23] . We will later refer to 
the theory behind the bang-bang version of UDD to derive always-on CAFE sequences. 

2.1. Dynamical Decoupling by Order in T 

Our task is to develop a class of control functions /3(f) so that the system evolution of 
the qubit over time interval f € [0, T], given by 

W(0,t) = Texp (~*J H(t )dr ^ , (5) 

closely approximates the identity operation on the qubit (to within a complex phase 
factor) at the final time f = T, where T is the time-ordering operator. We 
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approach this problem by expanding the interaction picture propagator in the Magnus 
expansion, as follows: 



where 



W int (0,T) = exp l-iJ2^k J, (6) 
V fc=o / 

fii = / n int {s)ds, (7) 
Jo 

^2 = 1 [ [ [n int ( Sl ),n int {s 2 )}d Sl ds 2 , (8) 

* Jo Jo 

and so forth. To simplify analysis, we transform the integration variables to 
dimensionless form as 

2 

Uj =- Sj -l, (9) 

so that the domain of integration is now [—1,1]. As a result, il n cx T n , so the expansion 
in ^ is a power series in T. 

We say that a sequence which satisfies il 3 ■ ~ for j = 1, . . . , n "decouples to 
n -order" because the first non-vanishing error is of order T n+1 ||W in t|| n+1 . We use 
approximate equality because the CAFE sequence we will derive is ultimately limited 
by numerical precision and physically realizable control fields, which may violate strict 
equality of the "decoupling constraints." However, we defend the introduction of any 
approximation errors with our simulation results. This approach to the dynamical 
decoupling problem is akin to time-dependent perturbation theory, and is useful when 
the noisy environment is not strongly perturbing the system and hence the Magnus 
expansion (in the interaction picture defined above) converges rapidly. 

The "zeroth-order" constraint which must be satisfied is that the control field 
itself (in absence of any bath effects) is identity. This could be modified, if desired, to 
give some single-qubit rotation, but we do not pursue this possibility here. Hence, 

P(T) = (mod27r). (10) 
This additional constraint must be appended because the initial set of constraints was 
derived in the interaction picture. For simplicity, we will use /3(T) = strictly in this 
work. 

The first- and second-order decoupling constraints can be derived from ^ and Q 
by nullifying any terms which are not identity operations on the qubit. For example, 
the first-order constraints (cx T) are 

J sm[(3(u)]du = 0, (11) 
-i 

cos[(3(u)]du = 0. (12) 
The second-order constraints (cx T 2 ) can be derived from the Magnus expansion as 
J J {sin[/3(iti)] - sm[(3{u 2 )}} d Ul du 2 = 0, (13) 

{cos[/?(ui)] - cos[/3(w 2 )]} duidu 2 = 0, (14) 

-l 

Ul 

sin [P{u x ) - (3(u 2 )] du x du2 =0. (15) 

-l 
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These continuous DD constraints are equivalent to those studied elsewhere, for 



example in Ref. 24 



2.2. A Continuous Control Function that Mimics Uhrig Dynamical Decoupling 

Finding arbitrary time-dependent control functions that satisfy the continuous DD 
constraints, even at first-order, is not trivial. We begin our search by studying the 
cumulative pulse area /3udd(£) for the Uhrig dynamical decoupling sequence (UDD) 
because this sequence is known to decouple by perturbative order using bang-bang 
pulses 29 . By transforming to the domain u £ [—1,1], the time-dependent rotation 



Ml 



in a UDD sequence of N pulses is given by 

[■u n 

/3udd(u) = 7T / ^5 



N 

i=i 



cos 



_JTT_ 

N + 1 



du\ 



J7T 

N + l 



(16) 



where S(u) is the Dirac delta function and Q(u) is the Heaviside step function. 

The integrated control function for UDD, /?udd(w), can be well approximated by 
a smooth function of the form 



&(«) = (JV+lJcOB-^-u), 



(17) 



where N corresponds to the number of pulses in the UDD sequence, which is plotted 
in figure fl] for N = 6. Note that Ref. 29 makes a convenient variable substitution, 



which is to operate in a transformed coordinate basis given by 9 = cos 1 (— u). In this 
basis, P s is proportional to 9. We utilize this transformation later. By differentiation 



of (171, one can determine that the corresponding control function is given by 

iV- 



a s (u) 



1 



2yT 



(18) 



Moreover, the interaction picture Hamiltonian can now be expressed in terms of 
Chebyshev polynomials when using this continuous control function since 

cos[(iV + 1) co8 -1 (-u)] = T N+1 (u), (19) 



sin[(iV + 1) cos -1 (-u)] = y/l-u 2 U N (u), 



(20) 



where T/v+i is the (N + 1) -order Chebyshev polynomial of the first kind and Un 
is the N th -ovdei Chebyshev polynomial of the second kind. The latter is particularly 
noteworthy because Vl — u 2 is the weight function for Un, so that \/l — u 2 Un(u) is 
orthogonal to any polynomial in u of order N — 1 or lower on the interval [—1, 1]. This 
property allows us to assert that the dynamical decoupling constraints containing just 
a "sine" function, as in (11) and (13), are made zero by this choice of function [29]. 



We can summarily state this as 



dui . . 



du m sm.[l3 s (uk)] = 



(21) 



for m < N — 1 and 1 < k < m [29] . Furthermore, half of the corresponding "cosine" 
constraints in the it-domain are zero by parity rules, since T/v+i( u ) nas definite 
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parity. The remaining constraints are not zero, so this control function has a nonzero 
constraint at first- or second-order. Figure [I] compares sin[/3 s (u)] to cos [/?udd(w)], 
also known as the UDD "switching function" y(t). This indicates that (3 S preserves 
many of the decoupling properties of UDD, as supported by numerical simulations to 
be discussed later. 

Before proceeding to derive new sequences based on smooth approximations 
of UDD, we note that the connection of UDD's switching function to Chebyshev 
polynomials allows an alternative interpretation of its decoupling abilities. To arrive at 
a switching function, we treat the time-dependent operator B(t) as a random classical 
function, following standard procedures surrounding the Born-Markov approximation. 
In a particular trajectory, a qubit in state co|0) + ci|l) evolves to co|0) + exp(i</>)|l) 
with the accrued phase during the UDD sequence 

B{t)y(t)dt. (22) 



o 



This integral may be rewritten 



2{N+1) ^(-iy Wj I PAruir. ,28) 



where the Tj are the locations of UDD pulses, i.e., Tj = Tsin 2 (j7r/[2(./V + 1)]), or, 
on the u axis, Uj = — cos(jx/(N + 1). The weights are given by Wj = ir/2(N + 1) 
for j = 0, N + 1 and Wj — w/ (N + 1) otherwise. These weight functions are those of 
Gauss-Lobatto-Chebyshev quadrature, which says that for order M, 



M 

E 

3=0 



^>rjin;) = j JM=du + R M , (24) 



where Rm is a remainder term which vanishes quickly with M and Uj = — cos(j7r/M). 
Hence we associate M, the order of Gauss-Lobatto-Chebyshev quadrature, with N + l, 
where N is the number of 7r pulses in a UDD sequence. The other term we must 
consider is (— This may be written as (— 1) J = T/v +1 (uj), evident from the 
formula T n (u) — cos(ncos _1 (w)). Therefore, if we define 

/(«) - {N + 1)T T N+1 (u) f U B[T(v + l)/2]dv, (25) 



then (23) is exactly in the same form as the left-hand-side of (24). We may therefore 
replace our sum with an integral with only higher-order corrections. Then, integration 
by parts and a return to the t domain yields 

, 4 




dtB(t)U N [-~l)J-[l--)+R N+1 . (26) 



Recalling that y/(t/T)(l - tJT) = \ZT — u 2 is the weight function for the orthogonal 
polynomials Un(—u), we see that UDD effectively extracts the iVth term of an 
expansion of B(t) in Chebyshev polynomials of the second kind. Lower-order 
components of this expansion are decoupled. Higher-order components still appear 
in the remainder term Rj^+i- 

This result provides two major conclusions. First, it provides an alternative 
viewpoint as to how UDD achieves decoupling: UDD eliminates lower order Chebyshev 



polynomial components of the bath field. Second, direct comparison of (26) 
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t/T 



Figure 1: (a) A smooth approximation to the UDD pulse area over time using a 
continuous control function. At prescribed instants in time, the UDD sequence applies 
a 7r-pulse, seen here as a jump in the cumulative pulse area. The smooth approximation 
intersects with the UDD sequence at these pulse times. The traces are divided by ir to 
show the jumps at each 7r-pulse as integral steps in this plot, (b) Comparison between 
the UDD switching function (blue) and the function sin(/3 s (u)) = vU— u 2 Un{u) 
for f3 s (red), where Un(u) is the A^ th Chebyshev polynomial of the second kind. 
Because yl — u 2 is the weight function for this polynomial, it has many of the same 
orthogonality properties as the UDD switching function. 



and p2| shows that, up to a remainder Rn+i which vanishes for high N, the 
discontinuous switching function y(t) is well approximated by the continuous function 
(4/ic)U N (2t/T - l)sJ(t/T){l-t/T) = (4/tt) sin [&(£)], as shown in figure^. This 
connection motivates our effort to find smooth versions of the UDD sequence. 

2.3. Variational Parameters for Improved continuous DD Control Functions 

We will now develop a method for improving the dynamical decoupling properties of 
the continuous-control sequence given by f3 B . We consider a technique for expanding 
the difference /3 B — /?udd as a Fourier series, where we subsequently truncate the 
series and tune the remaining terms for optimal performance. Recent work by Uhrig 
has adapted the UDD sequence to replace instantaneous pulses with time-dependent 



control functions, while preserving many of the error-suppression properties 25 
This followed a mathematical structure derived by Yang and Liu [29] used to prove 
that the original UDD sequence works as intended. Recall the variable substitution 
6 = cos _1 (— u). In the 0-domain, /3udd(#) is a "staircase" function while j3 s (0) — NO 
is a linear approximation to /?udd- The difference between the two sequences is a 
periodic sawtooth wave, which has a Fourier series given by 



1 

3=1 J 

The first attempt at improving the control function would be to truncate the 
Fourier series with the first several terms (as shown in figure [2]), but this results in 
only a modest improvement in decoupling. However, there is an insight to be gained 
here, because it turns out that any "sine" Fourier components given by X p sm(Np8) 
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Figure 2: Construction of the CAFE sequence. The left panel shows /3tjdd(*) in black, 
a staircase function with steps of tt. Using (28), the red curve closely approximates 



/3tjdd(£) as a sum °f A (i) and a 41-term Fourier series for a sawtooth wave (only 20 
terms of which are nonzero). The green curve truncates this series to 5 terms, only 
two of which are nonzero; an offset is added to the figure for clarity. The blue curve 
is the numerical solution for the 5-parameter CAFE^ 3 ' 5 ' sequence, found using the 
green curve as a starting guess, again with added offset. The right panel shows the 
corresponding a(t) functions for these approximations; the UDD sequence is a series 
of ^-functions shown as black lines. The continuous approximations, again offset for 
clarity, have poles at the endpoints. 



do not affect the "sine" DD constraints, in the sense that 



P 71 


rO m -i 


/ d9\ . . 


. 1 d6 m sin 


fa 


Jo 



= 



(28) 



for 1 < k < m and any real coefficients A p (see Appendix A). We can exploit this 
property to nullify some of the other DD constraints by varying the A p parameters 
without violating the above "sine" DD constraints. In particular, the "cosine" DD 
constraints do depend on the A p parameters, so we set up a system of nonlinear 
equations involving the first m constraints that are not already zero by identity and 
the first m variational parameters (A p ). These equations involve nested integrals over 
complicated functions, and so finding genuine roots to the equations can be numerically 
intensive. 

The optimized set of parameters we examine more closely is given by five DD 
constraints that are not identically zero for the control function 



5 

£ 

k=l 



Afc sin(4fc#), 



(29) 



where Ai . . . A5 are variational parameters we can change to satisfy the five DD 
constraints. As discussed previously this functional form automatically nullifies "sine" 



DD constraints as in (11) and (13). Furthermore, this form has definite parity (even), 
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Fourier 


Optimized Search 


Ai 





0.0017 


A 2 


1 


0.9121 


A 3 





-0.2869 


A 4 


y 2 


1.3520 


A 5 





0.4920 



Table 1: The first five A p parameters for a truncated Fourier series and the optimized 
sequence found by numerical solution of five decoupling constraints. 



so half of the "cosine" constraints can be eliminated by parity rules. The system of 
equations we solve for is given by: 

cos[/3(wi)]d«i = (30) 



dui 



du 2 sm[j3(u 1 ) - /3(u 2 )] = 



(i*i) 2 cos[/3(ui)]dui = 



-l 
l 



dui / du 2 (u 2 y 



sm 



(ui) - /3(u 2 )] = 



(u±) cos[/3(ui)]diti = 0. 



(31) 



(32) 



(33) 



(34) 



Importantly, the parity rules noted in Section |2.2| only apply to equations in the 
u-domain, not in rea l time. Still, this set of equations nullifies all of the 1 st - and 2 nd - 

. Moreover, several error terms at 3 rd - 



order constraints in ( 1 1 
are also canceled. 



15) 



4th 



and 5 -order 



The X p parameters were treated as variables to solve the system of equations in 



(30 - 34 1. The initial approximate solution was taken to be the truncated Fourier series, 



and Table [T] lists the initial parameters and their values after optimization. Figure [2] 
shows the optimized continuous control function with five X p parameters. We call 
such decoupling sequences with optimized parameters CAFE. To specify precisely 
which CAFE sequence is being considered, we denote each as CAFE^'™ 1 ) where N is 



the order of the continuous approximation to UDD, as appears in (17) and subsequent 



references to (3 S like (29). The quantity m is the number of variational X p parameters. 
For example, figure |^shows the control function a^(t) for CAFE^ 3 ' 5 ). 



2-4- Piecing Together Sequences: Splice and Invert 

The previous section described how to improve sequences by inserting m variational 
parameters and searching for a set that minimizes m previously nonzero DD 
constraints. The resulting control function oS m \t) has the undesirable property of 
having poles at the start and endpoints, which are impractical to implement in an 
experiment. However, the function is very steep here as well, so it seems reasonable 
to "chop off" this section of the sequence with minimal impact on the decoupling 
properties. This process is illustrated in figure [3a[ The problem here is that one must 
also ensure that the net action of the control function is the identity gate. Otherwise, 
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Figure 3: The CAFE^ 3,5 - 1 sequence can be made finite by removing the beginning and 
the end of the sequence, (a) The CAFE( 3 - 5 ' sequence "chopped" at the second zero 
from the end of the sequence, (b) Two iterations of the CAFE^ 3,5 ' 2 ) base sequence, 
showing that the original was chopped at the second and second-to-last zeros. The 
second time window is the negative of the first. This function is finite and smooth at 
all points. 



the sequence will decouple the qubit well, but the state of the qubit will be rotated 
deterministically around the Bloch sphere. This implies that the pulse area in f3(t) 
that is removed must be equal to 2ir. The pulse area can be calculated readily, but 
ensuring this area is precisely 2ir in a real experiment seems unreasonable, given how 
large the first derivative of a(t) is at the endpoints of the sequence. 

A rather simple and effective solution is to chop at zeros of a(t) near the start /end 



of the sequence, as shown in figure 3a For a(t) being an even function of time, the 
zeros lie symmetrically about the sequence, so we can splice together a window of 
time with control a(t) followed by an equal window of time with control — a(t) which 
has been chopped in the same manner; for a(t) being an odd function of time, one 
must time-reverse control in one of the windows. Figure |3b| shows how to use this 
method to make a simple bounded sequence by splicing and inverting two segments 
of the CAFE sequence from figure [2] The resulting sequence has several advantages: 

(1) it is finite and continuous, and the first derivative is continuous, at all points; 

(2) the total pulse area is automatically zero since the second half of the sequence 
reverses any action in the first half; (3) this sequence preserves much of its dynamical 
decoupling properties and performs quite well in simulation. We expand the notation 
for this decoupling function to CAFE^' 1 ™' 1 ") xL where r is the root-number at which 
the sequence is truncated (counting from the end) and L is the number of times the 
CAFE base sequence is repeated. For example, figure 3b shows a CAFE' 3 ' 5 ' 2 - 1 x2 
sequence. 
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3. Comparing different sequences 

We evaluate the performance of the CAFE sequences in comparison to a very similar 



set of sequences, namely the RUDD sequences of Uhrig and Pasini [24 25 . However, 
we make two alterations to the original RUDD sequence to make the comparison as 
appropriate as possible. First, the RUDD sequence as originally derived [24] begins 
and ends with very sharp pulses of area 2tt. Our numeric studies as well as those 



of Uhrig and Pasini 25 have shown that these pulses have very little effect on the 
decoupling abilities of the RUDD sequence. We therefore omit these pulses in all 
comparisons shown in this paper. Another modification we make is to invert the 
sign of each subsequent pulse in the RUDD sequence. We denote a RUDD sequence 
modified in this way by RUDD*-A-d, where N is the number of pulses and d is the 
duty cycle. 

The inversion of subsequent copies of CAFE was introduced to assure continuity, 
but this process also introduced a further advantage of adding robustness to low- 
frequency error in the amplitude of the control field. Such errors will be discussed 
in more detail in Sec. [6j In the case of RUDD with perfect pulses, the modification 
of inverting each pulse has no effect on the performance of the sequence. However, 
it introduces the same robustness to low-frequency control noise, and so therefore 
provides a closer comparison to the CAFE sequences. 

In this paper, we will limit our discussion to CAFE^ 3 ' 5,r ^ xL sequences, varying 
only the root number r and the number of repetitions L. We will compare these 
sequences to RUDD*-A-<i sequences of different pulse numbers N and duty cycles d. 
In general, increasing L or reducing r improves the decoupling ability of CAFE, while 
increasing N or reducing d improves the decoupling ability of RUDD*. Any of these 
actions, however, increases the maximum slew rate, max t \da/dt\, of the sequence. For 
RUDD* a high duty cycle d means that pulses are stretched to have longer duration, 
which reduces slew rate since integrated pulse area is constant (7r-pulses). Limits on 
slew rate are likely to be the primary reason for using continuous sequences such as 
CAFE and RUDD*, and therefore we use this criterion to choose which CAFE (3,5 ' r) xL 
sequence to compare to which RUDD*-A-g? sequence. In particular, for numeric 
studies we will look at two groups of sequences, each of roughly equal maximum 
slew rate. One group has a lower slew rate and is summarized in figure [4j a second 
with higher slew rate is summarized in figure [5| 

4. Decoupling a qubit from a quantum bath 

Our first evaluation of the CAFE sequence considers a fully quantum mechanical 
description of a bath. For this we employ the canonical central-spin problem, in 
which our qubit is coupled to 6 spins which are themselves coupled via dipole-dipole 
interactions, i.e. 

B Q = 4 J ^ ^ -g (35) 

B z = Y i C j Il (36) 

3 

where lj is the spin vector (using / = 1/2) for the bath spins and Yjk is the spatial 
vector connecting them. The cubic geometry we choose is shown in figure [6k. This 
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Figure 4: Four sequences with comparable maximum slew rate max ( T 2 \da/dt\ rs 
7 x 10 5 . The top row shows the continuous control field a(t) in units of T. The 
bottom row shows 0(t), the integral of a(t); this function would look like a square 
wave from to tt for an ideal 7r-pulse sequence such as UDD. 
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Figure 5: Three sequences with comparable maximum slew rate max t T 2 \da/dt\ ss 10 7 ; 
see figure [4] caption. 



type of simulation tests the DD sequences against a system with more degrees of 
freedom [c.f. (11-151] than the Born-Markov classical bath simulations we present 
later. However, these additional degrees of freedom also complicate analysis, so we 
present no analytic expectation for the result. We demonstrate only that CAFE 
successfully decouples a qubit from a fully-quantum environment in a manner highly 
comparable to RUDD*. 

The plot in figure [6]d shows the projection of the qubit Bloch vector in the 
interaction picture after averaging over the bath spins. Specifically, this simulation 
initializes into state 

Pi(0) = ^y^®^, (37) 

where a J is either a Pauli matrix or single-qubit identity (j = I, x, y, z) and 1 is the 
bath identity matrix. Numeric integration then solves dpj/dt = — i[H, pj] for each 
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Figure 6: Simulation of decoupling a qubit from spin lattice, (a) Diagram of the 
spin lattice. The qubit is surrounded by six bath spins in a cubic lattice configuration, 
(b) Qubit Bloch vector components as a function of time during decoupling sequences, 
in the interaction picture, for varying initial conditions. The four sequences plotted 
are those in figure [4j The grey curves correspond to no control. For this figure, 
a x is abbreviated X, etc. Note that in the case of no control, a qubit beginning in 
state p z (0) will not change, hence the lack of grey curve on the logarithmic plot of 
Tr[Zpz(t)]. This particular simulation used Cj = l/T for all bath spins j. 



j = x,y,z, and the three qubit directions k = x,y, z are plotted as Tr[cr Pj(t)] each 
time t for varying initial conditions j. 

To evaluate infidelity, we consider the procedure of quantum process tomography 
(QPT), in which the final state of the qubit, after the sequence is complete at time T, 
may be written 

P( T )= E XikoSpWo*. (38) 

j,k={I,x,y,z} 

In single-qubit QPT, measurements of the three components of the Bloch with four 
initial conditions — p x , p y , p z , and p- z — are sufficient to construct the matrix Xjk m 



general 30 , but since this system features unitary evolution prior to the final trace, 
projections for initial condition p_ z are exactly the negative of projections for initial 
condition p z . Further, for our particular interaction and bath, the Ising-like symmetry 
assures Tr[a v p x (t)} = Tr[a* p x (t)] = and Tr^p, (*)] = for all j ^ x. Hence the 
four choices of projection and initial condition shown in figure [6}d are sufficient to 
construct Xjk- The infidelity for dynamical decoupling is defined as 1 — Xnl which 
may be interpreted as the probability that something other than identity happened 
to the qubit. In general, for a single qubit in which the process is unitary evolution 
involving a bath followed by tracing that bath, one may readily show 

\U~ E Tr[aV,(T)]l. (39) 

j={x,y,z} 

Figure [7] shows the infidelity as a function of J, the overall strength of the dipole- 
dipole coupling within the bath, for both the low-slew-rate sequences and high-slew- 
rate CAFE and RUDD* sequences shown in figures [4] and [5] Varying J results in 
varying the correlation time of the bath. At low J, the bath is very slow, and here 



Xii 
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1/JT 1/JT 

Figure 7: Infidelity when decoupling from quantum bath with various sequences. The 
process infidelity is plotted as the overall timescale of the bath dipole-dipole coupling, 
1/J, is varied. (Left) Low slew-rate sequences shown in figure [4] (Right) High 
slew-rate sequences shown in figure [5] The points filled in with black are those 
corresponding to the trajectories of figure [6| 



dephasing without DD is maximal and determined entirely by the number of possible 
configurations of the bath spins. As J increases, even without DD the infidelity 
decreases due to motional narrowing, in which spin-flips within the bath reduce its 
coupling to the qubit. With CAFE and RUDD* sequences, it is clearly seen that 
DD strongly reduces the infidelity at low J, corresponding to long correlation times. 
The efficacy of decoupling reduces, at varying rates depending on the sequence, as J 
is increased; improved fidelity at even higher J is again due to motional narrowing 
effects. 

Figure [6] shows that RUDD* sequences and CAFE sequences are qualitatively 
different in how they decouple in time; in particular RUDD remains pulsed, while 
CAFE is fully continuous. However, figure [7] shows that they have roughly the same 
decoupling power per slew rate as JT is varied. To analyze their response to varying 
bath correlation times more quantitatively, we shift to a classical, Born-Markov type 
bath and a filter function analysis. 

5. Filter function analysis 

When studying the effects of noise on information, often it is instructive to 
analyze the system in the frequency domain. Dynamical decoupling schemes are 
effectively high-pass filters, meaning they suppress low-frequency components of the 
system-environment coupling while not reducing high-frequency noise. The "cut-off 
frequency" between these regions is often closely related to the rate at which one can 
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apply control pulses, but the precise relationship is less clear in examples like UDD, 
where pulses are not equidistant. To accurately describe how a decoupling sequence 
protects against noise with a given spectral density, we characterize each sequence 



by a "filter function" 31 , which characterizes the degree to which noise at a given 
frequency is suppressed. Filter functions are a meaningful way to compare different 
DD sequences because they provide insight into which method is best suited to a given 
noise spectrum. 

We require that the filter function be a linear response which is independent of 
the magnitude of the environment noise, but this is only appropriate under certain 
conditions. First, we invoke the Born-Markov approximation and represent the 
dcphasing bath with a scalar function B(t), so that the Hamiltonian is now 

H(t) = B(t)a z + a(t)a x , (40) 

where B(t) is a stationary random noise field with zero mean and spectral density 
S(oS) given by 



S(w)= dt(B(t)B(r)} cos(w(i - r)). (41) 

J — oo 

Second, we assume that the magnitude of the bath noise is sufficiently weak that we 



can calculate the qubit infidelity as a linear function of S(u)) (see Appendix B |. Under 
these conditions, we may then define the filter function F(z) by stating that the qubit 
infidelity X — 1 — Xiii or total error, is given by 

^S(.)F(.T). (42) 

The filter function for a given decoupling sequence can be calculated numerically 
in simulation, by modeling the noisy environment with a sinusoid. By measuring 
the degradation in qubit fidelity averaged over an ensemble of random phases for the 
sinusoid at frequency u), one obtains the linear response F(ljT). For control fields 
without a cr z -componcnt, the filter function may be theoretically expressed as 

2 



F(loT) 



ujT / e^ Tu+ ^ Tu ^du 



(i 



(43) 



For a general derivation of this expression in the continuous case, see |Appcndix B| 
This expression clearly reduces to the canonical expression for the filter function 
derived for sequences composed of ideal 7r-pulses such as UDD 23 s 31 . For such 
sequences, e 1 ^^ = cos[/3(i)] = y{t) is the switching function. Note that if the smooth 
approximation for y(t) corresponding to UDD, i.e. (A/Tv)U N (2t/T-l)y/(t/T)(l - t/T) 



as discussed in Sec. |2.2| is substituted into (431, the integral may be analytically 
performed to yield F(wT) « 16(iV+ 1) 2 J^ +1 (wT/2) for UDD-7V, an expression which 
yields a very close approximation to the exact UDD filter function, as stated without 
derivation in Ref. [23] , 

Figure [8] shows filter functions comparing different versions of CAFE and RUDD, 
calculated both with explicit integration and via Monte Carlo simulation. Examining 
figure |HJ we see that despite dissimilar time-domain construction, CAFE provides an 
intermediate limit between low-iV, low-duty-cycle RUDD sequences and high-iV, high- 
duty-cycle RUDD sequences. The sharp dips in the filter function are quite different 
between the two cycles, but these play little role in realistic noise bath models which 
typically integrate out such features. We turn to such models next. 
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Figure 8: Filter functions calculated from simulation (open circles) and numeric 



integration (solid lines) following ( 43 1 . The "noise" in the Born-Markov approximation 



is a sinusoid, so that sweeping frequency over many trials of the decoupling sequence 
reveals the degradation in fidelity, and hence the filter function F(u). 



6. Decoupling a qubit from Born-Markov classical baths 

To further analyze the effectiveness of CAFE sequences in comparison to RUDD* we 
simulate the total evolution of arbitrary qubit states coupled to classical dephasing 
baths. As in the filter function analysis, we invoke the Born-Markov approximation 



leading to the scalar noise field B(t) in (40) with noise spectral density S(lo) given 



by (41 1. We now consider bath noise with both Gaussian and Lorentzian spectral 
density. The reason we consider a classical bath is that by doing so, we can readily 
define parameters such as absolute magnitude and correlation time for a scalar noise 
spectral density, as opposed to the quantum bath which has far too many degrees of 
freedom to permit simplified analysis. 

Our simulations incorporate two sources of noise: environmental dephasing noise 
and noisy control fields. Control noise is inevitable in a practical experimental setting. 
Even a sophisticated decoupling sequence could yield mediocre results if such control 
noise is not taken into account; moreover, a DD sequence could possibly introduce 
more noise than it corrects. Due to these possibilities, we investigate the effect of 
control noise on CAFE and RUDD* by simulating their decoupling effectiveness with 
and without noise added to the control functions. 

Figure [9] shows the behavior of our example CAFE and RUDD* sequences under 
Born-Markov noise with a spectral density of Gaussian form, i.e. S(tu) cx exp[— (wt c ) 2 ], 
both with perfect control and with 1% control noise. Figure [To] shows a similar 
simulation with the spectral density of Lorentzian form, i.e. S(u>) cx 1/[1 + (ojt c ) 2 }. 




Figure 9: Final state fidelity for decoupling sequences, with environment dephasing 
noise treated as Born-Markov classical bath with Gaussian S(lu). Simulation results 
are shown with and without noisy controls, for both low slew-rate and high slew-rate 
sequences. The solid lines show numeric integrals of the numerically integrated filter 
functions (solid lines of figure |8| with the spectral density targeted by the simulated 
dephasing noise. No modification is made to the solid lines for the plots showing 
control noise. 
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Figure 10: Similar simulation results to figure [9| except with environmental dephasing 
noise treated as Born-Markov classical bath with Lorentzian S(ui). 



The control noise is modeled as a time-dependent fluctuation on the amplitude of the 
control field a(t), taken with root-mean-square magnitude 1% of a(t) and a Gaussian 
spectral density with correlation time T/2 for the plots shown. Theoretical curves, 
found by numeric integration, are shown for comparison; the approximations inherent 
in these curves fail most prominently at very low and very high correlation times r c . 
In the case of no control noise, CAFE's performance is roughly intermediate between 
the RUDD* sequences, as already evident from the filter functions analysis of figure [8] 
One should note the strong similarity between the results of figure [9] with no control 
noise and those of figure [7j the nearly identical behavior gives a clear example of a 
fully quantum bath behaving identically to a classical one, supporting the frequent 
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use of the Born-Markov approximation to describe spin baths. 

As control noise is added, we see that all sequences roughly follow their expected 
behavior at short correlation times, where decoupling is already least effective. In 
the regime of long correlation times, where DD performs best, a noise-floor is reached 
which depends on the sequences; this noise-floor is consistently better for the RUDD* 
sequences. It is tempting to associate this lower noise floor with a lower duty cycle for 
the RUDD* sequences with respect to CAFE's 100%, but this conclusion is questioned 
by the highly similar behavior of RUDD* sequences of very different duty cycle. As 
discussed in Sec. |3j the noise floor of the RUDD sequences without the modification 
of alternating the sign of pulses, not shown, is much worse. This simple method of 
alternating pulse signs has certainly added some robustness, but many improvements 
in optimization against control-noise are certainly possible in future work. 



7. Conclusion 



The CAFE sequence originated from the analysis of UDD as if it were continuous. 
Using principles from numerical quadrature, we showed that UDD's essential behavior 
could be captured by continuous approximations to its normally discontinuous 
description. The continuous approximation to UDD provided a heuristic for 
continuous sequence design. The RUDD sequences follow a similar path, but are 



based on a different heuristic 24|25 . Both sequence families can be extended to higher 
order. A promising feature of these types of sequences appears to be that, for these 
two families at least, sequences can be locally altered to eliminate short-time spiking 
behavior without appreciable detriment to their decoupling ability. Although different 
heuristics were employed, the two sequence families behave remarkably similarly, as 
we have demonstrated using a variety of analyses. 

The reason analysis of continuous-control decoupling is more complicated than the 
analysis of "bang-bang" sequences is because noise can disturb the qubit during control 
operations. More than a complication of analysis, however, this situation is a reality 
of experiment. We believe, based on our results, that filter functions remain relevant 
for continuous sequences, and they provide an intuitive picture for what components 
of the environment noise are suppressed. Furthermore, Monte Carlo simulations of 
semi-classical noise bath spectra and fully-quantum spin lattices indicate how well a 
sequence performs at protecting quantum information. By adding noise to the control 
field, we can test another source of error, the experimenter and the experimenter's 
equipment, and how well the sequence copes. We hope the heuristics underlying the 
construction of CAFE and the characterization methods we have used may be useful 
to future research in continuous dynamical decoupling. 
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Appendix A. Proof that Fourier series components do not affect 
"sine"-like DD constraints 

We begin with the definitions from before that 

ft(t) = (JV + lJcOB-^-tt), 

which in 0-domain (0 = cos _1 (— it)) becomes 

p s (e) = (n + i)6. 

We can now define a control function with m A-parameters given as 



(A.l) 
(A.2) 
(A.3) 



p a is a control sequence with m variational parameters designed to cancel m DD 
constraints not zero for (3 S alone. Here we show that this construction will not affect 
the "sine" -like DD constraints given by 



dti 



dt m sm[f3(t k )}=0 



(A.4) 



for 1 < k < m. Using integration by parts, the above set of constraints can be reduced 
to the system of equations 

rT 

t n sm(/3)dt = (A.5) 

for n = 0, 1, ... , (m — 1). Working in the 0-domain, this set of constraints is equivalent 
to 



f [1 - cos(0)] n sin(/3) sm(9)d6 = 0, 
Jo 

which is also equivalent to 

cos(p0) sin(/3) sin(0)d0 = 



for p = 0, 1, . . . , (m — 1). By substituting the expression in (A.3 1, 
cos(p0) sin[/3^ m) ] sin(0)ef0 = 



(A.6) 
(A.7) 
(A.8) 



cos(p9) sin(0) sin[(iV + 1)0] cos 



cos(p0) sin(0) cos[(iV + 1)0] sin 



A fc sin((JV + l)k6) 

k 

^A fe sin((7V + l)fc0) 



dO 



dQ. 



Note that cos[^ fe A& sin((AT + l)fc0)] can be written as the infinite sum 
Sfco a j C0S ((N + l)j0) for some real coefficients a,j since cos(0) is an even function 
in <p. Similarly, sin[^ fc A& sin((7V + l)fc0)] can be written as Y^jLi bj sin((A + l)j0) 
for some real coefficients bj since sin(</>) is an odd function in (j>. The DD constraints 
arc therefore satisfied if 



cos(p0) sin(0) sinL4(AT + 1)0] cos[B(iV + 1)0]<20 = 



(A.9) 



for integers A,B satisfying A > 0, B > 0, and p < N. This equation can be directly 
evaluated to always be true under these conditions. 
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Appendix B. Derivation of filter function for continuous sequences 

In this appendix we derive a general expression for the filter function of a continuous 
DD sequence. We generalize slightly from the main text, allowing a Hamiltonian of 
the form 

H{t) = 2eB(t)S z + A(t) ■ S. (B.l) 

Here A is a vector of control fields, one for each component, and S is the spin vector, 
i.e. 5 J = ai /2. The unitlcss parameter e is in place to keep track of orders of bath 
interactions, since the filter function is defined only at the lowest nonvanishing order. 
We use e to define that the filter function F(z) in terms of the infidelity X at vanishing 
e. Again, X is defined defined as 1— xn for the QPT-derived matrix due to evolution 
at finite e. Then F{z) is defined so that 

lim £ - 2 Z = J ^S{w)F{uT). (B.2) 

This type of e-based definition is not required for pulses composed of (5-function pulses, 
but in the continuous case where noise occurs during the control, the limit as e — > 
is critical for the linear-response theory upon which the filter function is based. 

We begin by transforming to an interaction picture, for which we need the control 
unitary Uo(t) — Texp[— i J Q A(t) • S]dr. In terms of this unitary, our interaction 
picture Hamiltonian is H(t) = U$(t)H(t)U (t). If A(t) chan ges in direction over 
time, Uo(t) in general may describe complicated trajectories over the Bloch sphere. 
In this most general case, it is not clear how to arrive at an analytic form for a filter 
function, and one generally must rely on numeric solvers. Here, we consider the simple 
case that A(t) = 2a(t)h, for constant direction described by unit vector n. The CAFE 
and RUDD sequences considered in this paper have all used n along the x direction; 
our definition of infidelity and assumption about the bath interaction indicate that any 
direction along the Bloch sphere equator is equivalent. In this appendix, we generalize 
slightly by allowing a constant z-component to n as well. Then, 
H{t) = 2eB(t){S z cos[/3(t)] + n z n • S{1 - cos[f3(t)}} - S ■ z x nsin^(t)]}, (B.3) 
using (3(t) as defined by 

We consider a single-qubit density matrix in the interaction picture 

PW = I+T.P'W, (B.4) 

3 

which evolves via H(t), to lowest nonvanishing order in time-dependent perturbation 
theory with ensemble averaging and standard Born-Markov separation of timescales, 
as 

' ,P =~ f 'dr{[H(t),[H(T),p(r)}}). (B.5) 



ii 



dt 

Correspondingly each Bloch-vector component pP evolves according to the following, 
where each repeated superscript is summed: 

dp> 



dt 



t dr2TT{SmH(t),[H(r),S k }})}p k (r) 



c 2 



r dLO 



o 

t 



= -e 2 



dr / — cos[u{t-T)]S{uj)M M {T)M ll {t)p k {T) (B: 
2ir 

,I T f ^cos[ W (^r)]5HP fe (i,r)/(r), 
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(B.7) 



M ke (t) = 2{e zkl cos[/8(i)] + n z h m e mM {l - cos[/3(f)]} 
+ [h k S ze -S kz h e ] sin[(3(t)}}. 

Here £ jW is the fully asymmetric Levi-Cevita tensor density and & k is the identity 
matrix. 

Now, we time-integrate p 3 to lowest order in e, to find 



pi(T) = 



dt / dr 



duj 
2^ 



cos[to(t-T)}S(cu)T jk (t,T) 



/(0) + O(e 4 ). (B.8) 



We use this expression in QPT, which in general gives 

1 " 



Xn (T) 



2 + 2p%(T)+2pl(T) + p z z {T) 
- PUT) - P X Z (T) pl(T) - p*_ z (T) p y_ z (T) 



(B.9) 



where p k = Ti{a k pj(T)} for initial condition p±j(0) — (a 1 ± oi)j2. As discussed in 
Sec.|4j the unitary evolution of this model assures pP z = —p : L z , simplifying the fidelity 
expression to 

l + pi(T)+py(T)+p z z (T) 



XI i (T) = 

As a result, the infidelity is simply written 

2 pT ft p 7 

I=- dt dr -cos[w(t-r)]%)Tr{r(f,r)} + 0(e 4 ), 
4 Jo Jo J 2tt 

and our generalized filter function is therefore simply 

f du [ dvcos[z(u- v)]Tr{T(Tu,Tv)}. 
Jo Jo 



F(z) 



(B.10) 
(B.H) 
(B.12) 



Upon taking the trace, the matrix r jfe (<, r) simplifies substantially to 

Tr[r(*,r)] = 8{cos[^) - /?(r)] + [n^] 2 [l - cos[/3(*) - (3(r)}}}. (B.13) 

In this paper, we only consider control fields for which fi z (t) — at all times. In this 
case our general expression for F(z) is 

F{z) = 2z 2 ( du ( dv cos[z(u - v)] cos[/3{Tu) - f3(Tv)}. (B.14) 
Jo Jo 

This easily simplifies to (431 in Sec.[5j 
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